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Abstract 

Remote observations of coronal holes have strongly implicated the 
resonant interactions of ion-cyclotron waves with ions as a principal 
mechanism for plasma heating and acceleration of the fast solar wind. 
In order to study these waves, a WKB (Wentzel-Kramers-Brillouin) 
linear perturbation analysis is used in the work frame of the collisionless 
multi-fluid model where we consider in addition to the protons a second 
ion component made of alpha particles. We consider a non-uniform 
background plasma describing a funnel region in the open coronal holes 
and we use the ray tracing Hamiltonian type equations to compute the 
ray path of the waves and the spatial variation of their properties. 

1 Introduction 



The ultraviolet spectroscopic observations made by SUMER (Solar Ultravi- 
olet Measurements of Emitted Radiation) and UVCS (Ultraviolet Corona- 
graph Spectrometer) aboard SOHO indicated that heavy ions in the coro- 
nal holes are very hot (hotter than protons by at least their mass ra- 
tio, i.e. Tj /Tp > mi /nip) with high temperature anisotropy (see, e.g.. 



Kohl et al . 1997 : Wilhelm et all 19981 ). This result is a strong indication 
for plasma heating through ion-cyclotron resonance (i.e., collisi onless en- 
ergy exchange between ions and wave fluctuations, see [Stb3 . 11992I . Chap 10) 
involving ion-cyclotron waves that are presu mably generated in the lower 
corona from small-scale reconnection events (lAxford and McKenzie , 1 19951) 
or locally throug h plasma micro-instabilities ([Mecheri and Marschl . 120071 : 
Markovskii I2OOII ) or turbulent cascade of lo w-frequency MHD-type waves 



towards high-frequency ion-cyc lotron waves (jLi et all Il999l : iHollwed . l200d : 



Ofman. Garv. and Vifia^ . 120021 '). This notion led to renewed interest in mod- 



els involving ion heating by high-frequency ion-cyclotron waves (llsenberg. Lee, and Hollwee 



,2000; H ollwee 



2000 



Marsch and Tul. 



2001 



Vocks and Marsch. 2001 



2004 : iBourouaine. Vocks. and Marschl . I2OO8I ) . For a detailed 



Xie. Ofman. and Vinasl . 



review on res- 



onant ion-cyclotron interactions in the corona, see iHollweg and Isenberg 
(I20p2|) and 1 1 plas ma properties and physical processes in coronal holes. 



see 



Cranmerl (|2009l '). Therefore, to to allow for a good physical interpreta- 



tion of the observed data, it became relevant to understand the properties 
of these waves in the naturally multi-ion and highly non-uniform coronal 
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plasma, particularly in preparation to the future data which will be pro- 
vided by the recently approved high resolution ESA-NASA mission "Solar 
Orbiter" (planned launch in 2014). 



Multi-ion linear mode analysis has been previously investigated bv lMann. Hackenberg. and Marsch 
( I997I ) in the frame work of the multi- fluid model but in the simple case of a 
uniform background plasma. Our aim in this paper is to go beyond this pre- 
vious study and consider a realistic background plasma with typical coronal 
profiles of density and temperature and a typical 2D open magnetic field 
model describing a funnel region in a coronal hole. We give a detailed study 
of ion-cyclotron waves by performing a Fourier plane wave analysis using the 
collisionless multi-fluid model where in addition to electrons (e) and protons 
(p), we consider a second population of ions, namely alpha particles (He^'*', 
indicated by a) with a typical coronal abundance. While neglecting the 
electron inertia, this model permits the consideration of ion-cyclotron wave 
effects that are absent from the one-fluid magnetohydro dynamics (MHD) 
model. Considering the WKB (Wentzel-Kramers-Brillouin) approximation, 
we first solve locally the dispersion relation and then perform a non-local 
wave analysis using the ray-tracing theory. This theory allows to compute 
the ray path of the waves as well as the spatial evolution of their properties 
while propagating through the funnel. 

This paper is structured as follows. In Sect. [2l we present the background 
plasma configuration we use in term of density, temperature and the 2-D 
analytical funnel model describing open magnetic field region in a coronal 
hole. Then in Sect. [3l we describe how the local and non-local (ray-tracing) 
linear mode analysis are carried out using the multi-fluid model. The results 
are presented and discussed in Sect. [H and finally our conclusions are given 
in Sect. [5l 



2 Background plasma configuration 



As a background plasma de i isity and temperature we use the model of 
Fontenla. Avrett. and Loeser ( 19931 ) for the chromosphere and Gabriel ( 19761 ) 
for the lower corona (see Fig. [H). The 2-D potential-field model derived 
by Hackenberg. Marsch. and MannI ( 20001 ) is used to define the background 
magnetic field (Fig. [T|) representing a funnel. Analytically, the two compo- 
nents of this model are: 
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Figure 1: Left: Electro nic density (wp.) and temperature (T ) model-profiles 
of the chromosphere ( Fontenla. Avrett. and Loeser . 19931 ) and the lower 
corona (|Gabriei[l97^ ^. Right: F unnel magnetic field geometry as obtain ed 
from a 2-D potential field model ( Hackenberg. Marsch. and Mann . 20001 ) . 



arctan- 



cosh^ sin^ + sin{ 



nd 
2L 



2ttx \ " 
L > 



sinh^ cos^ 



(1) 



with the following typical relevant parameters: L = 30 Mm, d = 0.34 Mm, 
Boo = 11.8 G, and B^ax = 1-5 kC. 

3 Basic equations 



The cold collisionless fluid equations for a particle species j are: 

^ + V-(n,v,) = 0, 
^j^ji-Q^ + • Vvj) - qjUjiB + Vj X B) = 0, 



(2) 
(3) 



where mj, nj, qj and Vj are respectively the mass, density, electric charge 
and velocity of a species j. Subscript j stands for electron e, proton p or 
alpha particle a (He^"*"). The electric field E and the magnetic field B are 
given by Faraday's law: 

VXE=-^. (4) 



3.1 Linearization procedure 

The linear perturbation analysis is performed by expressing all the quantities 
in the above equations as a sum of an unperturbed stationary part (with 
subscript 0) and a perturbed part (with subscript 1) much smaller than the 
stationary part: 

rij = noj{z) + nij,B = Bo(x, z) + Bi,E = Ei, Vj = 

with: nij <C noj, |Bi| ^ |Bo| . (5) 
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Considering the quasi-neutrality (i.e., J2j Qj'i^Oj = 0)) ambient electric 
field (i.e., Eq = 0) and no background velocity (i.e., vqj = 0). The zero- 
order terms cancel out when Eq. ([5]) is inserted into the Eqs. ©-(jH)- 
Neglecting the nonlinear products of the first-order terms, we get a system 
of coupled linear equations: 

iuj—^ — ik • vij = 0, (6) 

iuj^ij + % + vij X ^ j = 0, (7) 

ikxEi = ia;Bi, (8) 

where all the perturbed quantities have been expressed in the form of a 
Fourier plane wave. This Fourier analysis turns all derivatives into algebraic 
factors, i.e. d/dt — > —ioj and V — ?> ik, where a; is the wave frequency and k 
the wave vector. The quantity Vtj = qjEo/mj is the cyclotron frequency of 
species j. 



3.2 Dispersion relation 

To derive the dispersion relation, the above linearized equations have to be 
combined in order to obtain a linear relation between the current density Ji 
and the electric field Ei: 

3i = a-Bi, (9) 

where a is the conductivity tensor which is related to the dielectric tensor e 
through the following relation: 



e{u!, k, r) = I H (t(cj, k, r). 



(10) 



Finally, the loc al dispersio n relation is obtained using the theory of electro- 
dynamics (e.g.. IStixl . I1992I ): 



i:>(a;,k,r) = Det 



^k X (k X E) -Fe(a;,k,r) -E 



0, 



(11) 



where c is the speed of light in vacuum and r is the large-scale position 
vector. We choose the wave vector k to lie in the x — z plane, with k = 
k{sin6, 0, cosO). 



3.3 Principal waves properties 

In order to define certain wave properties we adopt the following coordinate 
system: e^^ and e^, are the unit vectors respectively in the direction of the 
ambient field Bq and the wave vector k (Fig. [TJ. We assume here that these 
two unit vectors lie in the x — z plane. The unit vector perpendicular to 
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efc, but still in the x — z, is e^p. The angle between Bq and k is (p. In this 
configuration, the phase velocity and the group velocity Vgr are given 
by: 

uj _ duj 1 d(jj 



Vph 



efc, 



A: '^^ ■ ^" dk"'' ' k 

The angle between the group velocity and the ambient field is given by: 



(12) 



if) = arctan 
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ok k am 



dk 
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(13) 



The helicity g (which represents the degree of circular polarization) and the 
electrostatic part of the wave ^ are given by (jKrauss-Varban. Omidi. and Quest! . 
I994I I: 

__,Re(^ e = -^. (14) 



1 + \Pa 



Pa 



where r epres ents the (/)-component of the polarization vector (see, e.g.. 



Shafranovl . ll967l ). 



3.4 Ray-tracing equations 

In the framework of the WKB approximation, the-ray tracing problem con- 
sists a t solving a system of ordinary differential equations of the Hamiltonian 
form ( Weinberg . 19621 ). These equations, which represent the equations of 



motion for the wave frequency w, the wave vector k, and the spa ce coor- 
dinate r, have been formulated by iBernstein and Friedland jigsi). In the 



simple case of a Hermitian dielectric tensor, they are given by: 

duj ^ dDiu;,k,r)/dt ^ 

dt dD{oj,-k,r)/doj ' ^ ' 

dk 3D(u;, k, r)/5r ^^^^ 



dt 5L>(w,k,r)/9w' 
dr 5Z)(u;, k, r)/9k 



(17) 



dt dD{io,k,r)/dco' 
Note that Eq. ()15p can be set to zero, because the dispersion relation does not 
explicitly depend on the time t (i.e., the background plasma is stationary). 
The above set of differential equations represents an initial-value problem 
which can be solved using initial conditions obtained from the local solutions 
of the dispersion relation Eq. ()lip . 



4 Numerical results 
4.1 Local wave analysis 

We consider a three- fluid cold plasma made of electron (e) , protons (p) and 
alpha particles (He^"*", indicated by a) with a density Ua = O.lup. The dis- 
persion relation Eq. ijlip is solved numerically at the funnel location (x=7.5 
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Two-fluid (e-p) Three-fluid (e-p-He^+) 




Figure 2: Top panel: Two- and three-fluid dispersion surfaces at the location 
(x=7.5 Mm, z=2.5 Mm) in the coronal funnel characterized by an inclination 
angle (p ~ 79° of the magnetic field with respect to the z-direction (see Fig. 
[1]). Here u) and k are normalized, respectively, to the proton cyclotron 
frequency and the inertial length VLp/VAp, with Vap representing the 
proton Alfven speed. Bottom panel: Two- and three-fluid dispersion curves 
at the same location and for different propagation angles with respect to 
the z-direction. 



6 




Figure 3: Basic wave properties at the location (x=7.5 Mm, z=2.5 Mm) 
in the coronal funnel characterized by an inclination angle (p ~ 79° of the 
magnetic field with respect to the z-direction (see Fig. [1]). The same nor- 
malization as in Fig. ([5]) is used. These results are given for the case of 
a propagation angle 6 = 50° with respect to the z-direction. The wave 
properties are : (a) refractive index kVAp/i^, (b) phase velocity Vph, (c) 
The group velocity Vgr, (d) helicity (degree of circular polarization) g, (f) 
electrostatic part of the wave ^, (e) angle between the direction of the group 
velocity and the ambient magnetic field ip. 
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Mm, z=2.5 Mm) which correspond to a hmit region between open and closed 
magnetic field lines where ion-cyc lotron could potentially be ge nerated from 

This loca- 
:i 79° with 



mi)- 



small scales reconnection events (jAxford and McKenzie 
tion is characterized by a magnetic field inclination angle of ip 
respect to the normal on the solar surface. For comparison purpose, the 
dispersion diagrams of a two-fluid (e-p) cold plasma are also presented on 
Fig.([2|). In this case, Eq.(|ll|) is a quadratic polynomial of degree 4, which 
means that 2 modes exist and each one is represented by an oppositely 
propagating (cj > and a; < 0) pair of waves. These modes, for which 
the dispersion curves and surfaces are shown on the left panels of Fig. ([1]), 
are the Ion-Cyclotron mode (IC) and the Fast mode. They are respectively 
the extensions of the usual Alfven and Fast MHD modes into the high fre- 
quency domain around w ~ fip where they are dispersive as compared to 
their non-dispersive MHD character at very low frequency <C ilp). In the 
case of the three- fluid model, Eq. ljlip is a quadratic polynomial of degree 6. 
Thus, 3 modes exist which in term of phase velocity ordering are: the Ion- 
Cyclotron mode 1, the Ion-Cyclotron mode 2 and the Fast mode. Indeed, 
the results represented on the right panels of Fig 12] reveals the appearance of 
a second ion-cyclotron mode IC2 in addition to ICl mode and the fast mode 
already present in the two-fluid model. For large k, the two IC modes reach 
a resonance regime at each one of the cyclotron frequencies, i.e. cj = Op 
and oj = ^la = rip/2. The heavy ion component strongly influences the 
dispersion branches at small frequency uj/Qp < 1 and small wave number 
kVyip/f2p < 1 (V^p is the proton Alfven speed). Indeed, we note the ap- 
pearance of a Fast mode cut-off frequency approxi mately at uJr.o/ ^p ~ 0.583 



in agreement with the algebraic relation given bv lMelrosel (|l986l ^: 



n 



Zq. 



1 + Aa{noa/nop) 

1 + Za{noa/nop) 



(18) 



with Aa = nia/mp = 4 and the charge number Za = 2. Also, at kVAp/^p ~0.5, 
the fast mode and the IC2 mode couple and convert at a frequency cor- 
responding to the so-called cross-over fr equency ujrr/^p ~ 0. 6 12, a lso in 
conformity which the relation derived bv lHollweg and Isenberg (|2002l ): 



nr. 



Za{na/np + Za/Al) 



1 ~l~ ZaUa/ flp 



1/2 



(19) 



On Fig. ([3]), we notice that the two IC modes starts at very low frequency 
w/fip <C 1 with a linear polarization (£» = 0) in accordance with MHD theory. 
For increasing k, they become very rapidly circularly polarized, right handed 
{q = +\) for the ICl mode and left handed [q = —1) for the IC2 mode. 
However, due to the coupling with the Fast mode, the ICl mode transforms 
rapidly into a left-handed circularly-polarized mode at large k. At large k. 
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these two modes reach a resonance regime at each one of the ion-cyclotron 
frequency uj = Qp and uj = J7q, = il.p/2, regime which is characterized by an 
infinite refractive indices kVyip/rip (see Fig. [3]). At very smah k, the IC2 
mode is characterized by a group velocity mainly parallel to the ambient 
field (i.e., ^ ~ 0) but slightly deviates from it | ^/j |~ 8° for large k (see Fig. 
[3]). Since the direction of the group velocity indicates the flow direction of 
the wave energy, this means that the IC2 mode transport its energy mainly 
along the background magnetic field. This mode is partially electrostatic 
with 0.23 < ^ < 0.33. On the other hand, at small k the ICl mode is quasi- 
electrostatic (i.e., ^ ~ 0) with an energy flow fairly oblique to the magnetic 
field with | '0 |~ 20° decreasing however to \ ip \^ 8° for large k. The Fast 
mode starts at the cut-off frequency cvco/^p = 0.583 at k = and is initially 
left-hand polarized {g = — !)• Owing to mode coupling with the ICl mode, 
it rapidly changes its sense of polarization around kYAp/^p = 0.5 to become 
a completely right-hand circularly-polarized wave {g = +1). The fast mode 
is weakly electrostatic ^ < 0.23 and its group velocity propagates rather 
obliquely to the ambient magnetic field (10° <[ ip \< 19°). 

4.2 Non-local wave analysis: Ray tracing 

In this section we go beyond the local treatment of the waves and perform 
a non-local wave study using the ray-tracing differential equations. They 
are solved employing initial conditions obtained from the local solutions of 
the dispersion relation Eq. ()lip presented in the previous section. The 
wave is launched at the initial location (xq = 7.5 Mm and zq = 2.5 Mm) 
with different initial angle of propagation, 6*0 = 20°, 40°, 60° and 80°. The 
computation have been performed for different initial values of the normal- 
ized wave number kg, but for the sake of consistency the results will only 
be given for ko=0.2. The conclusions which will be drawn are valid for all 
waves emitted at ko < 1. The ray trajectory of the waves while propagating 
in the funnel are shown on the bottom panels of Fig.Q and the spatial 
variation of their properties as a function of the z-coordinate are presented 
on Fig. ([5]) and ([6]). On the top panels of Fig.([H) are shown the polar plots 
of the group velocity for the three modes corresponding to different nor- 
malized wave numbers kVyip/r2p=0.2, 0.4,0.6 and 0.8. They reveal that the 
IC2 mode propagation is mainly anisotropic in the direction of the ambient 
field Bq and cannot propagate perpendicular to it. Thus the corresponding 
energy mainly flows along the magnetic field lines, with maximum group 
velocity at parallel propagation. On the other hand, the ICl mode and the 
fast mode propagate almost isotropically with no preferred direction, and 
consequently the energy flow is fairly isotropic. 

These results are confirmed by the ray-tracing computation which clearly 
shows that the IC2 mode is well guided along the field lines, since the ray 
paths (direction of the group velocity) for various initial angle of propagation 
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Figure 4: Top panels: Topology of the group velocity of the modes for differ- 
ent normalized wave number k and at the location x=7.5 Mm, z=2.5 Mm, 
with magnetic field Bq inclination-angle <p ~ 79°. The length of the radius 
at an angle of propagation 9 counted from the horizontal axis {9 = 0°) is 
equal to the group velocity. Bottom panels: The ray trajectory of the modes 
launched at the location (xo=7.5 Mm, zo=2.5 Mm) in the coronal funnel, 
with initial normalized wave number ko=0.2 and at different initial angles 
of propagation 9o = 20°, 40°, 60° and 80°. 
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Ion-Cyclotron mode 2 Ion-Cyclotron mode 1 



Fast mode 




Figure 5: Variation of the basic properties of the modes as a function of 
height z in the solar atmosphere starting at the launch location (xo=7.5 Mm, 
zo=2.5 Mm) in the coronal funnel (with a Bo-inclination angle <po ~ 79°) 
with an initial normalized wave number ko=0.2. These results are given 
for different initial angles of propagation ^o- The wave properties are : the 
refractive index kVAp/^, the phase velocity Yph and the group velocity V^^- 
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Figure 6: Variation of the basic properties of the modes as a function of 
height z in the solar atmosphere starting at the launch location (xo=7.5 Mm, 
zo=2.5 Mm) in the coronal funnel (with a Bo-inclination angle ipo ~ 79°) 
with an initial normalized wave number ko=0.2. These results are given 
for different initial angles of propagation ^o- The wave properties are : the 
helicity (degree of circular polarization) g, the angle between the direction 
of the group velocity and the ambient magnetic field -0 and the electrostatic 
part of the wave ^. 
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^0 nicely follow the magnetic field lines (bottom panels of FigS]). Indeed, as 
shown on Fig. ([6]), the maximal angular deviation is | "0 |~ 6.5° at the lower 
part of the funnel 3 < z < 5Mm whereas the propagation is quasi-parallel, 
i.e. tp 0°^ in its upper part. During this quasi-parallel propagation the 
IC2 mode is quasi-electrostatic, i.e. | ^ |~ 1, and is elliptically polarized 
with 0.1 <| |< 0.4. 

Contrary, the ICl mode and the Fast mode are unguided modes with their 
ray path having mainly a strait trajectory (bottom panels of FigU]). The an- 
gular deviation ^lJ between the direction of the group velocity and Bq, shown 
on Fig.©, varies between —50° and 70° for the ICl mode and —40° and 
60° for the Fast mode. The polarization is in general elliptical, right-handed 
(q > 0) for the ICl mode and left-handed {g < 0) for the Fast mode, except 
for 9 = 20° in the upper part of the funnel (z>= 6 Mm) where it is mainly 
circular right-handed {g = +1) for the ICl mode and circular left-handed 
= — 1) for the Fast mode. 

5 Conclusion 

Ion-cyclotron waves propagation in a coronal funnel have been studied via 
a linear mode analysis using the cold three- fluid (e-p-He^"*" ) model. 
First local solutions of the dispersion relation have been given for a defined 
region in the funnel. For comparison purpose dispersion curves in the case 
of a two-fluid (e-p) model have been also presented. They showed the pres- 
ence of two modes which are the extensions of the standard Alfven and Fast 
MHD modes into the high-frequency domain (w ~ Qp) where they are dis- 
persive. At large wave number, the Alfven mode (also named Ion-Cyclotron 
mode) experiences a resonance regime at oj = 0,p. In the three- fluid model, 
the results have shown that the consideration of a second ion population, 
namely alpha particles He^"*", strongly influences the dispersion branches as 
compared to the two-fluid model. Indeed, they are subject to mode coupling 
and conversion and show the appearance of a cut-off frequency (concerning 
the Fast mode). In addition to the first Ion-Cyclotron mode "ICl" and the 
Fast mode previously present in the two-fluid case, we note the appearance 
of a second Ion-Cyclotron mode "IC2" also subject to a resonance regime 
at w = at large k. 

The nonlocal wave analysis, which has been performed by solving the ray 
tracing equations, revealed that the IC2 mode is very well guided by the 
magnetic field since its ray path nicely follows the field lines independently 
of the initial angle of propagation. This is in agreement with the topology of 
the IC2 mode group velocity which predicted a fairly anisotropic direction 
of propagation in the direction of the magnetic field with a maximum value 
parallel to it. As a consequence, the energy associated with this mode can 
carried-out to the outer corona along the open magnetic field lines. On the 
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other hand, the ICl mode and the Fast mode are unguided since their ray 
paths do not follow the magnetic field lines. They can propagate in the 
funnel isotropically with approximately a strait trajectory. The basic wave 
properties in term of phase and group velocity, polarization and other im- 
portant quantities show in general a significant variation as the waves prop- 
agate in the funnel, revealing thus the importance of the non-uniformity of 
the medium. All these effects, and particularly the possible wave absorption 
and dissipation near the resonance frequencies of minor heavy ions might 
play a role in coronal heating. 
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